gr-qc/04 



Quantum resolution of black hole singularities 

Viqar Husainl and Oliver Winkler^ 

t Department of Mathematics and Statistics, University of New Brunswick, Fredericton, NB, Canada E3B 5A3 
* Perimeter Institute for Theoretical Physics, Waterloo, ON, Canada N2L 2W9 

(Dated: February 7, 2008) 

We study the classical and quantum theory of spherically symmetric spacetimes with scalar field 
coupling in general relativity. We utilise the canonical formalism of geometrodynamics adapted to 
the Painleve-Gullstrand coordinates, and present a new quantisation of the resulting field theory. 
We give an explicit construction of operators that capture curvature properties of the spacetime and 
use these to show that the black hole curvature singularity is avoided in the quantum theory. 
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It is well known that many solutions of classical general 
relativity have curvature singularities. The most com- 
monly encountered are the initial or final singularities in 
cosmological models, and the singularity inside a black 
hole. Just as the classical singularity of the Coulomb 
potential is "resolved" by quantum theory, it is believed 
that any candidate theory of quantum gravity must sup- 
ply a mechanism for curvature singularity resolution. 

This question has a long history with the first studies 
appearing in the late sixties using the canonical formula- 
tion of general relativity due to Arnowitt, Deser and Mis- 
ner (known as the ADM formalism) . The main approach 
has been (and continues to be) to first focus on specific 
models derived from general relativity by requiring that 
spacetimes have some prescribed symmetries. These re- 
ductions fall into two general classes, namely models that 
are particle mechanics systems, and those that are still 
field theories. The reduced model is then quantized in the 
standard Schrodinger approach, and the singularity res- 
olution question is studied by looking at operators that 
capture curvature information. If the spectra of such 
operators are bounded, it is taken as an indication of 
singularity avoidance. 

The earliest works on singularity avoidance concerned 
cosmological models. Inparticular the Gowdy cosmolo- 
gies received attention Q, Q, Surprisingly however, 
the black hole singularity problem has not been stud- 
ied to date in ADM geometrodynamics, inspite of rela- 
tively early formulations of the classical canonical the- 
ory in spherical symmetry 0, Q , and more recent work 
Furthermore, the spherically symmetric gravity- 
scalar field system is of central importance for a full study 
of Hawking radiation. A complete quantization of this 
system could give a solution to all the questions and puz- 
zles that arise from the semiclassical Hawking radiation 
calculation, of which singularity resolution is but one. 

Our focus in this paper is the asymptotically flat, 
spherically symmetric gravity-scalar field system. We 
first formulate the classical canonical field theory adapted 
to the Painleve-Gullstrand coordinates. This formu- 
lation appears the most relevant to the gravity-scalar 
field quantization problem, since there is no "extended" 



Schwarzschild solution with various sectors to worry 
about, and the coordinates are not singular at the hori- 
zon. This classical part of our work complements the 
Schwarzschild parametrization given in 0, 0] • 

For quantization, we use a non-standard approach 
based on recent work by the authors on singularity reso- 
lution in cosmology which in turn is inspired by the 
loop space representation , its application to cosmol- 
ogy jlflj - an d by work of Halvorson|llj. In this quan- 
tization the configuration and translation operators are 
represented on a Hilbert space, but the momentum op- 
erator (the infinitesimal generator of translations) is not 
defined. We construct a class of operators corresponding 
to phase space functions that are classically singular, and 
show that these operators are bounded. 

We begin by presenting an outline of the classical 
Hamiltonian theory of the spherically symmetric gravity- 
scalar field model. Since we would like to adapt this sys- 
tem to the Painleve-Gullstrand (PG) coordinates, let us 
recall that the black hole in these coordinates is given by 
the metric 
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The spatial metric e a b given by the constant t slices is 
flat and the extrinsic curvature of the slices is 
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where n a = x a /r for Cartesian coordinates x a . The 
black hole mass information is contained only in the ex- 
trinsic curvature, which in the canonical ADM variables 
{lab, K ab ) determines the momenta Tr ab conjugate to the 
spatial metric q a b- 

The phase space of the model is defined by prescrib- 
ing a form of the gravitational phase space variables q a b 
and TT ab , together with falloff conditions in r for these 
variables, and for the lapse and shift functions iV and 
N a , such that the ADM 3+1 action for general relativity 
minimally coupled to a massless scalar field 
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is well defined. The constraints arising from varying the 
lapse and shift are 

H = -J= (z ab Z ab - \^ ^qR(q) 

_ 47rG ( 1 P| + JqqrtdMfi = (4) 



C a = D c it c a -P <j> d a <i> = Q, 



(5) 



where 7r = ir ab q ab . The corresponding conditions for 
the matter fields (j) and P^ are determined by the con- 
straint equations. The falloff conditions imposed on the 
phase space variables are motivated by the PG metric 
JU , which itself is to be a solution in the prescribed class 
of spacetimes. These conditions give the following falloff 
for the gravitational phase space variables (for e > 0) 



q a b — e a b + 



r 3/2+e 

-i _ 9 a \0^) , h ab {6, 



0{r- 2 ) (6) 
+ 0(r" 2 ), (7) 
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where / ab , g ab , h ab are symmetric tensors, n ab — fr ab / Jq, 
and q = detg a &. 

In this general setting we introduce now the following 
parametrization for the 3-metric and conjugate momen- 
tum for a reduction to spherical symmetry: 

R(r,t) 2 

q a b = A(r, t) 2 n a n b H — (e ab - n a n b ) (8) 
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Substituting these into the 3+1 ADM action © shows 
that the pairs (P, Pr) and (A, P A ) are canonically conju- 
gate variables. We note for example the Poisson bracket 



Rf dr, e iXPR ^ \ = i2G\f(ry XPRir \ (10) 



which (below) will be the bracket represented in the 
quantum theory. The falloff conditions induced on these 
variables from J7J are the following 

R = r + Oir- 1 ' 2 -*), P R = r- 1 ' 2 + 0{r- 1 -') 

A = l + 0(r- 3 / 2 - e ), P A = r 1 / 2 + G(r- £ ) 

= r- 1 / 2 + 0{r~V 2 ~ e ), P$ = r x ' 2 + 0(r~ e ). 

These together with the falloff conditions on the lapse 
and shift functions 

N r = r~ 1/2 + 0{r- 1/2 - e ) N=l + 0(r- e ) (11) 

ensure that the reduced ADM 1+1 action 

Sr = dtdr (P R R + P a A + P ) 

+ / dtN r AP A (12) 



obtained by substituting © into © is well defined. This 
completes the definition of the classical theory. The two 
reduced constraints, the Hamiltonian and radial diffeo- 
morphism generators are first class. The surface term is 
needed to ensure functional differentiability of Sr. 

At this stage we perform a time gauge fixing using 
the condition A = 1 motivated by PG coordinates. This 
is second class with the Hamiltonian constraint, which 
therefore must be imposed strongly and solved for the 
conjugate momentum P\. This gauge fixing therefore 
eliminates the dynamical pair (A,P\), and leads to a 
system describing the dynamics of the variables (R,Pr) 
and (<fi, Ptj,) jl2l| . The reduced radial diffeomorphism gen- 
erator 



P A (P,^P R! P ) + 



PrR' = o 



(13) 



remains as the only first class constraint, which also gen- 
erates the local dynamics. 

To address the singularity avoidance issue, we first ex- 
tend the manifold on which the fields R etc. live to in- 
clude the point r — (ie. from R + to Rq ), which in 
the gauge fixed theory is the classical singularity. We 
then ask what classical phase space observables capture 
curvature information. For homogeneous cosmological 
models, a natural choice is the inverse scale factor a(t). 
This observable can be represented as an operator in both 
the connection-triad variable based loop representation 
[Io| , as well as in an ADM variable quantization [j| . For 
our model, a guide is provided by the gauge fixed theory 
without matter where it is evident that it is the extrin- 
sic curvature @ that diverges at r = 0, which is the 
Schwarzschild singularity. This suggests, in analogy with 
the inverse scale factor, that we consider the field variable 
1/P as a measure of curvature. A more natural choice 
would be a scalar constructed from the phase space vari- 
ables by contraction of tensors. A simple possibility is 



- 1 ( Pv Pr 
^2\R^ + KR 



(14) 



The small r behaviour of the phase space variables en- 
sures that any divergence in tt is due to the 1/P factor. 
We therefore focus on this. A first observation is that 
the configuration variables R{r, t) and <p(r, t) defined at 
a single point do not have well defined operator realiza- 
tions. Therefore we are forced to consider phase space 
functions integrated over (at least a part of) space. A 
functional such as 



Rf = 



drfR 



(15) 



for a test function / provides a measure of sphere size in 
our parametrization of the metric. We are interested in 
the reciprocal of this for a measure of curvature. Since 
P ~ r asymptotically, the functions / must have the 
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falloff f(r) ~ r ~ for Rf to be well defined. Using this, 
it is straightforward to see that 1/Rf diverges classically 
for small spheres: we can choose / > of the form / ~ 1 
for r << 1, which for large r falls asymptotically to zero. 
Then Rf ~ r 2 and diverges classically for small 

spheres. 

A question for the quantum theory is whether 1/Rf 
can be represented densely on a Hilbert space as 1/Rf- 
This is possible only if the chosen representation is such 
that Rf does not have a zero eigenvalue. If it does, we 
must represent 1/Rf as an operator more indirectly, us- 
ing another classically equivalent function. Examples of 
such functions are provided by Poisson bracket identities 
such as 



1 



R 



f\ 



iGf(r) 



-iP R {r) 



iP R {r) 



(16) 



as the functions / do not have zeroes. The representa- 
tion for the quantum theory described below is such that 
the operator corresponding to Rf has a zero eigenvalue. 
Therefore we represent 1 /Rf using the r.h.s. of (|16f) . The 
central question for singularity resolution is whether the 
corresponding operator is densely defined and bounded 

We now proceed with a quantization of the 1+1 di- 
mensional scalar theory of the fields R and 4>. A guide- 
line is the idea to turn an expression like e~ tp {x^ , e lp }, 
which resembles Eqn. I|16|) . into a quantum operator that 
has eigenvectors that correspond to the classical singu- 
larity. It turns out that the implementation of this idea 
requires that plane waves be normalizable states in the 
quantum theory. This renders a standard Schrodinger 
type quantization unsuitable. Instead one is forced to 
use a quantization based on an enlarged quantum con- 
figuration space. This type of quantization has already 
been fruitfully applied to an investigation of the initial 
big bang singularity in the FRW cosmological model 8] , 
where there is a discussion of some technical aspects of 
the quantization. We describe here the ingredients nec- 
essary for our discussion of the black hole singularity. 

As a prelude, let us consider the one-dimensional non- 
relativistic particle with Hamiltonian H(x,p). The basis 
states are |a;) = \exp(ixp)) with inner product (x\y) = 
5 Xj y, with the Kronecker delta indicating normalizability 
(as opposed to the case of a delta function in standard 
Schrodinger quantum mechanics). The enlarged quan- 
tum configuration space is the so-called Bohr compact- 
ification of K. (named for mathematician Harald Bohr; 
see |nj for details). Position operators act multi- 
plicatively, x\x) = x\x), while momentum operators exist 
only in their exponentiated form as finite shift operators 
e l vp\x) = \x — y). It should be stressed that the apparent 
discreteness evident here is not the same as the one due 
to introduction of a fixed space lattice at the classical 
level. 



The quantization of free scalar fields in this framework 
is similar, and has been developed in 0,0]. A quantum 
field is characterized by its excitations at N points in 
space. The important difference from standard quantum 
field theory is that here such states are normalisable. A 
typical basis state is 



■ a>N t ; °1 



■ bN 2 



(17) 

where the factors of L in the exponents reflect the length 
dimensions of the respective field variables, and Ofc, 6; are 
real numbers which represent the excitations of the scalar 
quantum fields R and <p at the radial locations {xk} and 
{yi}. The inner product on this basis is 

(ai . . . a Nl ; h, . . . bpf 2 \a[ . . . a' Nl ; b[ . . . b' N2 ) 



N 2 



if the states contain the same number of sampled points, 
and is zero otherwise. 

The action of the basic operators are given by 



Rf \ai... a Nl ;bi . . . b N2 ) = 
L 2 2_j a kf{xk)\a\ ■ ■ ■ a Nl ; b x 

k 



e iA i P Ji (xj)| ai _ _ ,a Nl ;bi 

t&l . . . , CLj — Xj 



>N 2 



(18) 



(19) 



where aj is if the point Xj is not part of the original basis 
state. In this case the action creates a new excitation at 
the point Xj with value — Xj. These definitions give the 
commutator 



R he r\P R (x) = -\f( x )L 2 e lXP ^ x ). 



(20) 



Comparing this with l|10|) , and using the Poisson bracket 
commutator correspondence, we see that L — y/2lp, 
where I p is the Planck length. There are similar operator 
definitions for the canonical pair (</>, Pa)- 

Using the expressions for the basic field operators, we 
can construct an operator corresponding to a classical 
singularity indicator: 



\Rt\ 



|oi ...a Nl \bi 



_ e -iP R (xj) 



°N 2 ) 



Ipf(Xj) 



Rf\ 



?iP R (xj) 



x|<zi . . . a Nl ;h . . . b N2 ). 



(21) 



One can show that the basis states are eigenvectors of this 
operator, and that all resulting eigenvalues are bounded. 
We illustrate this with the state 



Ja P R (r=0)\ 



(22) 
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which represents an excitation ao of the quantum field 
Rf at the point of the classical singularity: 

R f \S ao ) = (2l 2 P )f(Q)a |S O0 ), (23) 

which is clearly bounded. This shows that the singu- 
larity is resolved at the quantum level. In particular if 
there is no excitation of Rf at the classical singularity, ie. 
a.Q = 0, the upper bound on the eigenvalue of the inverse 
operator is 2/1%. This fact also gives an interesting result 
concerning areas of the symmetry 2-spheres: The inverse 
area has an upperbound proportional to l/lp. 

Having established that the curvature spectrum is 
bounded at the kinematical level, let us ask if this conclu- 
sion is affected by dynamics. This involves considering 
the action of the remaining constraint (|13[1 on basis states 
(|17(l . Our strategy is to show that this action is well de- 
fined even on states of zero eigenvalue of Rf, which is the 
indication that quantum evolution continues without di- 
vergence to the location of the classical singularity. This 
property is also the content of "dynamical singularity res- 
olution" in earlier works on cosmology [a, 1HJ| . It shows 
that local quantum evolution is well defined, unlike its 
classical counterpart. 

The reduced constraint (|13(l and the evolution equa- 
tions derived from it ^| contain factors of 1/R, which 
are the only sources of classical divergence at R = 0. 
This constraint is quantised by representing these fac- 
tors by powers of the 1/Rf operator defined above. The 
other terms in the constraint contain the field variables 
and their conjugate momenta. The former are diagonal, 
and the latter are represented as shift operators [lq. In 
both cases their action on a basis state is finite. Thus 
the only source of classical divergence in the constraint 
gets resolved. This may be compared to the resolution 
of the Coulomb divergence in atomic physics where the 
expectation value of V(r) ~ l/r is finite and bounded 
above, although the details of the mechanism are differ- 
ent. There it provides a lower bound on energy. In our 
case it provides dynamical curvature singularity resolu- 
tion via finite action of the constraint operator on states 
with zero eigenvalue of Rf. 

The second aspect of dynamical singularity resolution 
concerns the surface term, which may be viewed as gen- 
crating dynamics at spatial infinity. This "true" Hamil- 
tonian is given by N r Pa(R,4>, Pr, P^) evaluated at the 
point r — oo (since the sphere integrals just give a 47r 
factor). As before, this contains a factor 1/R which is 
bounded above at the quantum level for the same rea- 
sons given above. Therefore (xp\P\\tj^) is bounded for 
basis states \ip), and so is the expectation value in time 

evolved states (if>\P\l/ 'RfP\\ijj). 

Let us summarise these results and compare them with 
what has been achieved for cosmology. We presented a 



quantization of the gravity-scalar field system in which 
there are operators corresponding to classical curvature. 
All such operators contain a 1/Rf factor, which we have 
shown to be bounded. We have shown that the upper 
bound on curvature is unaffected by time evolution via 
the remaining constraint, which encapsulates the quan- 
tum dynamics of the system. This establishes both kine- 
matical and dynamical resolution of black hole singular- 
ities. 

The work on singularity resolution in homogeneous and 
isotropic cosmology in both the connection [l(j and met- 
ric variable [8j approaches does not utilise a time gauge 
fixing. The action of the Hamiltonian constraint is stud- 
ied on a kinematical Hilbert space, and dynamical singu- 
larity resolution is established by noting that the quan- 
tum constraint equation allows evolution past the state 
corresponding to the classical singularity. What we have 
shown here is that the same holds true in a field theory 
setting with radial inhomogeneity, in a fixed time gauge. 
This requires careful definition of the relevant field the- 
ory operators, which is much closer to what is expected 
for full quantum gravity. 

Conceptually, our work rests on the use of PG coor- 
dinates in a field theory setting with matter. The use 
of PG coordinates permits a unified treatment of the 
singularity resolution problem in a matter coupled set- 
ting. In contrast quantising the source free black hole in 
Schwarzschild coordinates requires a separation of the in- 
terior " cosmological" and "exterior" asymptotically flat 
regimes 01 • This advantage will be useful for a study of 
the full quantum collapse problem. 

The quantisation also provides a setting for stud ying 
other operators such as null geodesic expansions [lfij . 
These permits a definition of black holes in quantum 
theory without recourse to classical "horizon boundary 
conditions" at the classical level, which has the potential 
to address Hawking radiation in a fully quantised setting. 
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